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INTRODUCTION 


The  work  presented  herein  was  motivated  by  a  collection  of  unusual  strain 
data  acquired  during  a  series  of  test  firings  of  a  120-mm  gun  tube  during  the 
latter  part  of  1985.  During  these  tests,  circumferential  strains  exceeding 
three  times  those  predicted  by  the  customary  Lam6  design  formula  (ref  1)  were 
recorded  from  the  outer  surface  of  this  tube  a  few  feet  from  the  muzzle. 
Subsequent  inspection  of  the  strain  data  showed  that  these  strains  increased 
dramatically  with  projectile  velocity  and  possessed  an  oscillatory  waveform  with 
a  frequency  approaching  15  kHz. 

It  was  soon  discovered  that  the  extremely  large  strains  were  caused  by  a 
projectile  velocity  which  was  very  close  to  critical.  The  existence  of  such  a 
critical  velocity  had  been  predicted  in  the  literature  at  least  thirty  years  ago 
(refs  2-4),  but  to  the  best  of  our  knowledge,  no  observations  such  as  these  have 
been  reported  in  the  literature.  In  the  following,  the  essence  of  the  'critical 
velocity'  theory  of  axisymmetric  vibrations  of  a  circular  cylinder  is  reviewed 
and  the  theory  is  extended  to  include  the  case  of  non-axially  symmetric  vibra¬ 
tions  which  can  be  excited  through  coupling. 

CRITICAL  VELOCITY  THEORY 

This  theory  predicts  a  limit  as  to  how  fast  the  tube  deformation — in  the 
immediate  vicinity  of  the  projectile — can  be  made  to  travel  before  some  sort  of 
wave  develops.  An  exaggerated  view  of  this  deformation  when  the  projectile 
velocity  is  low  is  shown  in  Figure  1.  Under  these  'quasi-static'  conditions, 
the  deviation  from  the  Lam6-predicted  deformation  is  less  than  3  percent. 
However,  as  the  projectile  velocity  approaches  a  certain  critical  value,  the 
tube  deformation  in  the  neighborhood  of  the  base  of  the  projectile — that  is,  at 

References  are  listed  at  the  end  of  this  report. 
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the  front  of  the  moving  pressure — grows  dramatically.  As  with  classical  reso¬ 
nance,  the  deformation  is  theoretically  unbounded  when  the  critical  velocity  is 
reached  unless  damping  is  present. 


DEFORMATION  OF  BORE  SURFACE 
(STATIC) 

Figure  1 

The  simplest  equation  containing  the  essential  physics  of  the  situation  is 
(ref  2) 

0  +  R*  w  +  ffl  B  =  QU-H<x-vt))  (1) 

where  Q  is  a  constant  and  represents  the  magnitude  of  the  pressure  which  is 

assumed  to  be  moving  at  constant  velocity  v.  H  is  the  heaviside  step  function 

H(x-vt)  =0  x  <  vt 
=1  x  >  vt 

In  this  equation,  w  is  the  radial  displacement  of  a  point  on  the  median  surface 
of  a  cylindrical  shell  located  at  a  distance  x  along,  and  R  from,  the  central 
axis;  h  is  the  shell  thickness  and  is  assumed  to  be  small  compared  to  R;  m  =  ph 
where  p  is  the  mass  density  of  the  shell  material;  D  =  Eh3/12(l-i>2) ;  E  is 
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Young's  modulus  of  elasticity;  v  is  Poisson's  ratio;  and  v  is  the  velocity  of 
the  moving  pressure,  assumed  to  be  finite  and  constant. 

Conventionally,  steady-state  solutions  to  Eq.  (1)  are  sought  under  the  con¬ 
ditions  that  the  displacement  remain  bounded  at  x  =  ±  «  and  that  the  stresses 
and  displacements  be  continuous  at  the  location  of  the  pressure  front  x  =>  vt. 
Usually,  a  change  of  variable  from  x  to  {  is  made  where  £  =  x-vt.  Then  the 
moving  pressure  front  (the  projectile)  is  always  at  5  =  0  and  the  partial  dif¬ 
ferential  equation  (Eq.  (1))  becomes  an  ordinary  differential  equation  which  is 
easily  solved.  In  particular,  these  solutions  have  the  form  w  =  Ae±-ik£  and  are 
steady  when  seen  by  an  observer  moving  with  the  pressure  front  at  (  *  0.  k  is 
the  wave  number  and,  in  general,  is  complex.  Only  when  k  is  real  does  the 
assumed  form  of  the  solution  represent  a  wave.  To  find  what  waves  can  exist 
naturally  in  the  cylinder,  one  sets  Q  =  0,  and  substitutes  the  assumed  solution 
into  Eq.  (1).  It  is  seen  that  real  waves  are  possible  for  those  values  of  k 
which  are  the  real  roots  of  the  equation 

v(k)  =  k /  ?  (1 +(£)*)  (2) 

where  y*  -  and  v  is  the  phase  velocity  (real).  A  plot  of  these  wave 

numbers  versus  phase  velocity  is  called  a  dispersion  curve  and  is  shown  in  Fig¬ 
ure  2  (Cc  *  KE/p).  This  plot  shows  that  waves  with  low  wave  number  (long  waves) 
travel  with  phase  velocities  which  decrease  with  wave  number,  while  those  with 
high  wave  number  (short  waves)  have  phase  velocities  which  increase  with  wave 
number.  This  happens  because  of  two  competing  restoring  forces  in  Eq.  (1).  The 
tube  can  deform  as  a  cylindrical  membrane  in  which  case  the  second  term  of  Eq. 
(1)  dominates  the  behavior,  or  the  tube  wall  can  undergo  axisymmetric  flexure  in 
which  case  the  first  term  dominates.  The  fact  that  these  two  restoring  forces 
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compete  to  produce  a  minimum  in  the  dispersion  curve  of  Figure  2  is  the  impor¬ 
tant  part  of  the  critical  velocity  theory. 


k 

Figure  2.  Phase  and  group  velocities — thin  shell  theory. 

For  the  purpose  of  discussion,  it  can  be  considered  axiomatic  that  if  waves 
are  generated  by  a  moving  axi symmetric  pressure  Q,  the  phase  velocity  of  these 
waves  will  be  the  velocity  of  this  moving  pressure.  That  is,  the  load  must  be 
in  phase  with  the  wave(s)  it  creates.  With  this  in  mind.  Figure  2  shows  that 
such  waves  are  possible  provided  the  load  velocity  equals  or  exceeds  the  minimum 
possible  value.  Let  us  assume  for  the  moment  that  the  load  velocity  (the  pro¬ 
jectile  velocityj  somewhat  exceeds  this  minimum.  According  to  Figure  2,  there 
are  two  waves  with  this  phase  velocity.  Now,  physically  we  know  that  the  energy 
contained  in  these  two  waves  must  radiate  away,  not  toward,  the  source  of  the 
disturbance,  namely  the  moving  pressure  front.  Further,  it  is  known  that  energy 
travels  not  at  the  phase  velocity,  but  at  another  velocity  called  the  'group 
velocity,'  related  to  the  phase  velocity  as  follows: 
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Vg<k)  =  k  ♦  v(k)  (3) 

where  v(k)  is  the  phase  velocity  (dispersion  curve)  from  Figure  2.  A  plot  of 
group  velocity  versus  wave  number  is  also  shown  in  Figure  2  and  provides  a  cri¬ 
terion  for  deciding  whether  a  given  wave  belongs  in  the  region  (  >  0  or  {  <  0. 
For  example,  the  wave  corresponding  to  the  point  p  has  a  group  velocity  which  is 
smaller  than  its  phase  velocity.  Since  the  wave  energy  must  flow  away,  not 
toward,  the  pressure  front  at  {  *  0,  this  wave  belongs  in  the  region  5  <  0  only, 
i.e.,  it  is  a  'trailing  wave.'  Conversely,  the  wave  corresponding  to  the  point  q 
is  a  'head  wave,'  belonging  to  the  region  ?  >  0.  Assigning  the  two  waves  to 
their  appropriate  regions  results  in  the  following  solution  for  the  midwall 
displacement  w(x,t)  when  the  load  velocity  v  is  greater  than  the  minimum  phase 


velocity  of  Figure  2 


-b 2 

b*~-~a 


;  cos  af  +  1  ;  £  <  0 


.(2)  -a* 

'  »  B*"-~a*  cos  b< 


;  5  >  0 


where 


1-t-i  .  / i.z.l  .  „  .  /l-i-6  .  /*.: 

o  r  o  w  r  o  r  o 


-  1  _  mv2  /r2 

2  :  =  T  EhD 


and  C  =  QR2/Eh  approximates  the  Lam6  displacement.  Note  that  for  \  »  1,  the 
solution  for  x  <  vt  approaches  2C,  twice  the  Lam6  displacement. 

It  can  easily  be  seen  from  Eq.  (3)  that  should  dv/dk  ever  vanish,  the  group 
velocity  and  the  phase  velocity  would  be  equal  and  energy  could  not  radiate  away 
from  the  pressure  front,  but  would  continually  build  the  deformation  in  the 
neighborhood  of  the  front  as  time  progressed,  i.e.,  resonance  would  result. 

Thus,  the  minimum  phase  velocity  of  Figure  2  is  indeed  a  'critical  velocity'  and 
it  is  the  near  attainment  of  this  velocity  which  caused  the  high  strains  in  the 


120-mm  gun  tube. 
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Finally,  if  the  velocity  of  the  moving  pressure  is  less  than  the  minimum 
possible  for  wave  formation,  the  wave  number  k  is  complex  and  the  solution  to 


Eq.  (1)  has  the  form  of  a  damped  harmonic 
w(l)  ed£  h*  _  c* 

_  *  -—  (_cos  C£  +  ----- —  sin  c?)  +  1  ;  ?  <  0 

(5) 

w(2)  <j*  _  c* 

_  =  (cos  c?  +  — 2cd"  sin  ;  5  <  0 


where 

and 

A  REFINED  MODEL  FOR  THE  120-MM  QUN  TUBE 

From  the  standpoint  of  gun  tube  design,  it  is  important  to  be  able  to  pre¬ 
dict  critical  velocities  as  accurately  as  possible  and  to  be  able  to  predict  the 
steady-state  deformation  at  any  velocity  of  the  pressure  front.  Thus,  it  is 
necessary  to  use  a  model  which  is  not  restricted  to  cylinders  with  thin  walls. 
Axisymmetric  equations  of  motion  for  thick-walled  cylinders  have  been  derived  by 
Mirsky  and  Herrmann  (ref  5)  and  are  considerably  more  complicated  than  Eq.  (1). 
They  are  used  to  obtain  the  results  which  follow  in  much  the  same  way  as 
discussed  previously.  These  equations  and  the  details  leading  to  their  solution 
have  been  reported  by  the  author  (ref  6). 

Although  transient  effects,  boundary  reflections,  non-uniformity  of  wall 
thickness,  and  variable  pressure  (projectile)  velocity  are  ignored,  steady-state 
calculations  for  thick-walled  cylinders  nevertheless  produce  results  in 
remarkable  agreement  with  measured  values  when  the  projectile  velocity  is  close 
to  critical.  (The  assumption  of  constant  projectile  velocity  is  justified  in 
the  forward  regions  of  many  gun  tubes  where  the  projectile  velocity/travel  curve 
is  relatively  flat.)  Figures  3  and  4  compare  predicted  and  measured  values  of 
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dynamic  amplification 


circumferential  strain.  Figure  3  shows  the  tube  deformation  in  the  neighborhood 
of  the  projectile  and  Figure  4  shows  the  maximum  strain  as  a  function  of  projec¬ 
tile  velocity.  The  strains  have  been  normalized  with  respect  to  the  static 
values  predicted  by  the  Lam6  formula.  The  disagreement  between  the  measured  and 
predicted  strains  at  low  velocities  may  be  partly  due  to  transient  deformation 
and  non-axially  symmetric  deformation  induced  through  coupling.  Neither  of 
these  is  accounted  for  in  the  preceding  axisymmetric  steady-state  theory.  The 
latter  is  the  main  subject  of  this  report. 

AXISYMMETRIC  AMO  NON-AXIALLY  SYMMETRIC  WAVE  COUPLING 

In  addition  to  the  axisymmetric  equations  mentioned  in  the  last  section, 
Mirsky  and  Herrmann  (ref  7)  have  also  derived  a  set  of  five  equations  governing 
both  axisymmetric  and  non-axially  symmetric  motions  of  a  uniform  hollow 
cylinder.  Gazis  (ref  8)  has  gone  a  step  further,  accomplishing  the  same  using 
the  more  accurate  three-dimensional  equations  of  elasticity.  These  investiga¬ 
tors  have  computed  the  dispersion  curves  for  several  of  these  modes.  Viewed  in 
the  context  of  critical  velocity  theory,  it  is  interesting  that  the  dispersion 
curves  of  the  first  non-axially  symmetric  mode  and  the  previously  considered 
axisymmetric  mode  have  minimum  values  in  close  proximity  to  one  another.  This 
is  shown  in  Figure  5.  Both  curves  were  constructed  using  the  equations  of 
Mirsky  and  Herrmann.  The  ordinate  values  have  been  normalized  with  respect  to 
the  velocity  of  shear  waves  in  the  material,  Cs  =  VG/p,  where  G  is  the  shear 
modulus  of  the  material.  The  abscissa  is  the  non-dimensional  wave  number  h/L, 
where  L  is  the  wave  length.  The  non-axially  symmetric  mode  referred  to  in  the 
figure  corresponds  to  a  beamlike  deformation  of  the  axis  of  the  cylinder.  Since 
the  moving  pressure  is  perfectly  axisymmetric,  the  only  means  for  exciting  real 
waves  of  this  type,  in  par  icular  the  resonance  indicated  by  the  minimum,  is  by 
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some  sort  of  coupling  whereby  energy  can  be  exchanged  between  the  axisymmetric 
and  non-axially  symmetric  modes. 


Figure  5.  Uncoupled  dispersion  curves  (e  =  0) — Mirsky  and  Herrmann, 
thick-wall  shell  theory. 

One  source  of  modal  coupling  which  is  of  importance  in  gun  tube  design  is 
the  non-uniformity  of  the  thickness  of  the  tube  wall.  Intuitively,  it  is  clear 
that  when  a  tube  having  an  eccentric  wall  thickness  is  dynamically  pressurized, 
the  bore  centerline  tends  to  move  as  well  as  the  tube  wall.  In  effect,  the 
eccentricity  of  the  wall  thickness  creates  a  non-axially  symmetric  wall  stiff¬ 
ness  and  inertia.  A  good  indication  of  what  may  be  expected  has  been  obtained 
by  neglecting  the  asymmetric  stiffness  and  including  only  the  eccentric  inertia 
coupling.  This  is  the  work  reported  herein. 

Basically,  it  is  assumed  that  we  are  dealing  with  a  cylinder  which  has  a 
sinusoidal  distribution  of  mass  around  its  circumference  and  that  this  distribu¬ 
tion  is  uniform  along  the  direction  of  the  bore  axis.  Thus 

P  3  P0(l  +  e  cos  9)  (6) 

represents  the  mass  distribution  of  the  cylinder  with  e  serving  as  the  eccen¬ 
tricity  parameter.  pQ  is  the  mass  density  of  the  material.  The  coupling  effect 
sought  does  not  require  the  use  of  sophisticated  thick-wall  cylinder  equations 
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such  as  those  used  by  Mirsky  and  Herrmann  or  those  employed  by  Gazis,  and  con¬ 
sidering  the  idealism  reflected  in  neglecting  asymmetric  stiffness  effects  and 
the  assumed  mass  distribution  of  Eq.  (6)  above,  the  use  of  such  equations  is 
not  required.  However,  the  essential  character  of  the  dispersion  curves — 
especially  near  the  dual  minimums  and  throughout  the  long  wave  regions — must  be 
contained  in  whatever  equations  are  chosen.  It  is  surprising  that  the  Donne! 1- 
Mushtari  equations  (ref  9)  do  not  suffice  for  this  purpose,  failing  to  provide  a 
dispersion  curve  for  the  beamlike  vibration  which  passes  through  the  origin. 
Consequently,  the  equations  of  Flugge  were  chosen.  Though  more  complicated  than 
those  of  Donnell -Mushtari ,  these  equations  provide  the  necessary  features  in  the 
dispersion  curves  and  are  considerably  simpler  than  those  of  Mirsky  and 
Herrmann. 

The  Flugge  equations  of  motion  are  as  follows  (ref  9): 
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{i-v)  a* 

o 

a9  ii-v]  a9 

u 

o 

~2~  30* 

3s9  ~2~  3s30* 

o 

3(i-v2  a* 

{3-v2  39 

v 

o 

2~  as* 

~2~  is*a0 

a9  ii-v)  a9 

p -v|  a9 

!  +  ,  §!. 

w 

o 

as9  ~2~  aia§* 

~2~  as*30 

1  30* 

y 

(7) 


where  k  =  h*/12R*,  s  =  x/R,  and  0  is  measured  tangentially,  u,  v,  and  w  are  the 
displacement  components  in  the  x,  0,  and  R  directions,  respectively.  Q  = 
pR* (l-v* )/Eh  is  the  applied  load  due  to  the  pressure  p. 

Solutions  to  Eq.  (7)  are  assumed  in  the  form  of  a  linear  superposition  of 
axisymmetric  motion  (u0,w0)  and  motion  from  the  first  non-axially  symmetric  mode 
(uj.vj.wj).  (vQ  can  be  omitted  because  it  produces  an  equation  for  torsional 
motion  which  is  uncoupled  from  the  non-axially  symmetric  variables.)  Thus, 

U(S,0,t)  =  U0 ( X , t )  +  U} (x, t )cos  0 
v(s,0,t)  =  vjfx.tjsin  0 

w(s,0,t)  a  WQ(x,t)  +  Wj (x, t )COS  0  (8) 

Equations  (8)  can  be  considered  as  approximation  functions  for  use  with  the 
Galerkin  method  (ref  10).  The  variations  <5u  =  6u0  +  du^  cos  0,  6v  =  5vj  sin  0, 
and  5w  =  5w0  +  6wj  cos  0  multiply  the  first,  second,  and  third  of  Eqs.  (8), 
respectively.  The  resulting  expressions  and  Eq.  (6)  for  p(0)  are  substituted  in 
Eq.  (7)  followed  by  an  integration  between  the  limits  0  =  0  and  2 tt.  This 
results  in  ‘He  five  equations: 
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33w0  3w0  e  3*u^  3*u0  3*u0 

is3”  +  V  3s”  ’  2  it*”  '  3t*"  +  is*" 


3*wi  3wj  k  3wj  iy+li  dvj  _  3*ua  3*uj 

is3"  3s"  “  2  U  1  3s  ?  3s  3t*  3s* 

,  (1-U)  (l-p)  3*uo 

-  k  ut  -  i---i  U1  -  e  *  0 

k  /*>  3*wl  a*vl  3k(l-w)  3*vl  (l-p)  3*vl  (P+1)  3ul 

2  is*”  '  W1  "  it*'  +  ii*'  +  "2““  ii*”  V1  2'"  di” 


3*w^  3*wc 


34wr 


3  4  Ur 


3ur 


2  it**  +  it*”  +  k  "«S  *  kw0  +  w0  -  k  55,*  +  v  g5"  s  Q 


*y.o  +  k  .  2k  ?i!!i  +  Wl  +  e  a!”o .  (3:p> 

!■*  K  3s4  3s*  W1  e  3t*  *  2  3s*  V* 


3t* 


3s* 


_  k  +  „  ?“2  .  k  11:^1  -.2  =  o 

3S3  3S  2  3S 


(9) 


Steady-state  solutions  to  Eqs.  (9)  can  be  obtained  in  exactly  the  same  way  as 
the  steady-state  solution  to  Eq.  (1).  Under  the  substitution  £  *  s  -  vt  (where 
t  =  [E/p0R* (1-p* ) ]^t) ,  Eqs.  (9)  become  five  ordinary  differential  equations  with 
£  as  the  independent  variable. 

•ii  .  €Va  > 1  11 

-kwQ  +  PWq  -  -y  uj  +  (l-v*)u0  =  0 

III  I  t  II  II 

-2kwj  +  [p(k+2)  -  k]wj  (1+p)vj  +  2(l-v*)uj  -  2ev*u0 
+  [(k+l)(p-l)]u1  =  0 

k(3-p)wj  -  2w^  -  [2v*  +  (3k+l)(v-l)]vi  -  2v^  -  (1+p)uj  =  0 


ev*  11  1 1 1 1  11  111  1 

-y  wj  +  kwQ  +  v*wQ  +  (k+l)w0  -  kuQ  +  pu0  =  Q 

1111  11  ' 1  k  11  5  •  * 

kw}  +  (v*-2k)wi  +  WJ  +  ev*w0  +  ^  (p-3)vj  +  v^  -  kuj 

+  (w-i)  +  y]ui  =  0  (io) 
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The  complementary  solutions  of  Eqs.  (10)  are  obtained  by  substituting  solutions 
of  the  form  Aje"ia^  for  each  of  the  variables  u0,  uj,  v^,  w0,  and  wj,  where  Aj 
is  a  different  arbitrary  constant  in  each  substitution  and  o  is  the  wave  number 
The  result  is  a  set  of  five  linear  equations  which  only  have  non-trivia!  solu¬ 
tions  if  the  determinant  of  the  coefficients  of  the  Aj  vanishes,  i.e., 


■  — 

a* (v2-l) 

a2ev* 

2 

0 

— 

-ia(v+a*k) 

r  ■  1 1  —  ~ 

0 

a2ev2 

a2 (v2-l)  +  \  (k+1) 

-  \  (k+1) 

-  y  (1+w) 

0 

-io(|-  +v+a2k 

-  -) 

2 ' 

0 

y  (1+u) 

|-  [2v2+v(3k+l ) 

-3k-l]-l 

0 

-y  (v-3)  -  1 

-ia(v+a2k) 

0 

0 

-a2 (v2-a2k) 

+k+l 

Hi 

0 

-ia  (~  +v+a2k  -  \) 

-y  (3-v)+l 

-a2ev2 

-a2 (v2-a2k-2k) 

+  1 

(11) 


Setting  this  determinant  to  zero  results  in  the  dispersion  relation.  The  first 
two  branches  are  shown  in  Figure  6  for  the  value  e  =  0.1,  The  similarity  to 
those  of  Mirsky  and  Herrmann  (Figure  5)  is  evident. 

The  dispersion  curves  of  Figure  6  represent  the  natural  (unforced)  waves  of 
the  cylinder  which  are  the  complementary  solutions  to  Eq.  (10).  As  was  the  case 
in  dealing  with  Eq.  (1),  the  stresses  and  displacements  are  required  to  be  con¬ 
tinuous  at  the  location  of  the  pressure  front.  In  terms  of  the  variables 
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defined  in  Eqs.  (8),  the  following  quantities  must  be  continuous  at  {  «  0  (ref 
11): 

h*  '  '  '  1 

4R*  (W1  +  vl>  +  V1  -  ui  , 

I  t  t 

V(Wi+Vi)  +  Uj  -  Wj 

f  I  I 

wO  “  uo  • 

v-1  '  ( u-3 )  ’  •  •  •  «  » 

( ~2_ ) u 1  +  (y-2)W!  +  --y-  VJ  +  Wj  -  Uj 

til  It 

«  Wq  -  Ug  f 

and 

■  i 

Vl.  «0<  wl>  *0'  W1 


Figure  6.  Coupled  dispersion  curves  (e  =  0.1)  using  shell  theory  of  Flugge. 

The  continuity  relations  plus  the  boundedness  criteria  at  infinity  and/or 
the  group  velocity  argument  are  sufficient  to  uniquely  determine  the  displace¬ 
ments  once  the  particular  solutions  have  been  determined.  The  particular  solu¬ 
tions  corresponding  to  the  load  terms  on  the  right-hand  side  of  Eqs.  (9)  may  be 
found  by  eliminating  all  but  one  of  the  dependent  variables  through  the  use  of 


Cramer's  rule  and  the  commutativity  of  the  differential  operators  (ref  12). 

This  procedure  is  much  simpler  than  inverting  Fourier  transforms — a  common  solu¬ 
tion  technique  used  in  connection  with  these  problems  (cf.  (ref  3)). 

From  the  previous  discussion  in  connection  with  Eq.  (1),  several  load 
velocities  are  of  interest.  For  example,  from  the  dispersion  curves  of  Figure 
6,  it  is  clear  that  a  real  wave  will  be  excited  regardless  of  how  small  the  pro¬ 
jectile  velocity  might  be.  As  a  typical  value,  v  *  0.382  Cs  has  been  chosen  and 
the  corresponding  mid-wall  radial  displacement  at  the  top  of  the  cylinder  (6  = 

0)  is  shown  in  Figure  7  using  the  value  e  =  0.1  (10  percent  mass  eccentricity). 

1 . 4 

O  1,2 

I-  1 

<r 

*  .8 

h- 

Z  .6 

LU 
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<r  .  2 

o 

» — « 

Q  -.2 

-.4 

-10  -5  0  5  10  15  20  25 

S-UT 

Figure  7.  Radial  displacement  from  coupled  Flugge  equations, 
v/Cs  =  0.382,  e  =  0.1. 

The  ordinate  in  Figures  7  through  10  is  the  ratio  of  the  dynamic  displacement  in 
the  radial  direction  (w)  to  the  static  axisymmetric  radial  displacement  as 
approximated  by  the  particular  solution  for  wQ.  This  radial  displacement  is 
essentially  quasi-static  except  for  the  existence  of  a  wave  which,  from  the 
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group  velocity  argument,  exists  in  front  of  the  pressure  only.  The  presence  of 
any  'head  wave'  will  be  felt  at  the  gun  muzzle  prior  to  the  arrival  of  the  pro¬ 
jectile  and  therefore  constitutes  a  possible  cause  of  round  inaccuracy  at  the 
target. 

Another  load  velocity  of  interest  is  v  =  0.515  Cs  which  is  just  below  the 
lowest  critical  value.  Again,  e  =  0.1  and  the  radial  displacement  at  the  top  of 
the  cylinder  is  shown  in  Figure  8a.  The  head  wave  of  Figure  7  is  evident.  The 
axisymmetric  portion  of  this  displacement  is  shown  in  Figure  8b,  so  that  the 
difference  between  the  two  figures  is  the  contribution  of  the  beamlike  motion  at 
this  velocity  and  is  substantial. 


-10  -5  0  5  10  15  20  25 


S-UT 

Figure  8a.  Radial  displacement,  9  =  0,  v/Cs  =  0.515,  e  =  0.1. 
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Figure  8b.  Axisymmetric  component  of  Figure  8a. 

The  presence  of  a  maximum  value  in  the  long  wavelength  region  of  the  lower 
dispersion  curve  is  also  cause  for  interest.  It  is  known,  however,  that  the 
dispersion  curve  for  the  next  non-axially  symmetric  mode  (in  which  the  tube 
distorts  more  or  less  into  an  elliptical  shape)  passes  through  this  region  and 
has  a  minimum  value  there.  This  mode  is  not  included  in  the  present  analysis 
and  therefore  the  model  is  considered  to  be  inaccurate  in  this  region. 

Probably  the  most  interesting  load  velocity  is  that  which  falls  between  the 
minimums  (vcri,vcr2)  of  the  two  curves  of  Figure  6.  For  a  load  velocity  exactly 
halfway  between  these  minimums,  i.e.. 


vcrl  +  vcr2 
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Figure  9  shows  that  the  axisymmetric  and  beamlike  motions  are  nearly  equal  in 
magnitude  for  values  of  e  as  small  as  0.01.  This  is  apparently  due  to  the  fact 
that  the  minimums  (critical  velocities)  of  the  two  branches  of  Figure  6  approach 
each  other  as  e  diminishes.  Consequently,  a  load  velocity  halfway  between  them 
is  closer  to  both  critical  values  when  e  is  smaller.  This  is  shown  in  Figure 
10 — a  blowup  of  Figure  6  in  the  vicinity  of  the  minimums  in  the  neighborhood  of 
the  critical  velocities  for  two  values  of  e.  From  one  point  of  view,  the  non- 
axially  symmetric  wave  acts  as  a  vibration  absorber  for  the  axisymmetric  motion. 
This  is  not  likely  to  reduce  the  maximum  stress  in  the  cylinder,  however,  since 
the  wall  will  then  be  subjected  to  biaxial  stresses  of  comparable  magnitudes. 
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Figure  9a.  Axisymmetric  component  of  radial  displacement,  9=0, 


v  = 


vcrl  +  vcr2 


=  0.541,  e  =  0.01, 
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